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The rheology of suspensions of Brownian, or colloidal, particles (diameter d < 1 /im) differs 
markedly from that of larger grains (d > 50 /im). Each of these two regimes has been separately 
studied, but the flow of suspensions with intermediate particle sizes (1 /im < d < 50 pm), which occur 
ubiquitously in applications, remains poorly understood. By measuring the rheology of suspensions 
of hard spheres with a wide range of sizes, we show experimentally that shear thickening drives 
the transition from colloidal to granular flow across the intermediate size regime. This insight 
makes possible a unified description of the (non-inertial) rheology of hard spheres over the full size 
spectrum. Moreover, we are able to test a new theory of friction-induced shear thickening, showing 
that our data can be well fitted using expressions derived from it. 


Complex fluids, polymers, colloids and surfactant so¬ 
lutions find wide applications, partly because of their 
highly tuneable behavior under deformation and in flow. 
The success of the mean-field ‘tube’ model for poly¬ 
mers jjLj, which describes how each chain is constrained by 
thousands of neighbours, means it has long been possible 
to predict ab initio their linear and non-linear rheology 
from the molecular topology with very few free param¬ 
eters. In particular, a scaling description is available of 
the dependence of rheology on molecular weight. 

However, progress in suspension rheology has been 
more difficult [5]. The small number of nearest neigh¬ 
bours (order 10) rules our any mean-field description: 
local details matter. It is now possible to predict the low- 
shear viscosity of a suspension of Brownian hard spheres 
(HS, diameter d < 1 pm) up to volume fractions of </> < 0.6, 
and the rheology of granular HS (d > 50 pm) is increas¬ 
ingly being studied. Surprisingly, however, how the rhe¬ 
ology of HS changes over the whole size spectrum re¬ 
mains unknown, because the behavior in the industrially- 
ubiquitous intermediate size regime, 1 < d < 50 pm, has 
not been systematically explored. We offer such an explo¬ 
ration in this Letter, and show that the physics bridging 
the colloidal and the granular regimes is shear thickening. 

The rheology of colloidal HS is well known [3H5] : the 
viscosity is determined by the particle volume fraction, 
if >, and the dimensionless shear rate, or Peclet number, 
Pe (= tb 7, the shear rate 7 non-dimensionalised by the 
Brownian time, tb, needed for a free particle to diffuse 
its own radius). At Pe « 1 the flow is Newtonian; the 
viscosity becomes immeasurably large at </> g » 0.58 01. 
Shear thinning starts at Pe < 1, reaching a second New¬ 
tonian regime at Pe » 1 with a viscosity that diverges 
at random close packing [2], 0 rcp ~ 0.64, the densest 
amorphous packing for lubricated (frictionless) HS. 

Since tb scales as d 3 , granular HS inhabit the Pe » 1 
regime at all practical shear rates. Extrapolating naively 
from the above description of colloidal flow, one expects 
Newtonian behaviour with a viscosity diverging at <^rcp • 
Experiments do find a Newtonian viscosity, but it di¬ 


verges at a cf> that is lower than 0 rcp [7], the precise 
value being dependent on interparticle friction [S], 

How suspension rheology transitions from colloidal to 
granular behaviour as d increases has not been theo¬ 
retically predicted, and remains experimentally unclear. 
Previous experiments on size dependence mm stayed 
in the colloidal regime, or used highly polydisperse sys¬ 
tems near the granular limit nni- Thus, no unified de¬ 
scription over all particle sizes is yet available. 

We study the rheology of the intermediate size regime 
using poly-methylmetlracrylate (PMMA) spheres steri- 
cally stabilised by poly-12-hydroxystearic acid (PHSA) 
with d = 268 to 3770 nm, dispersed in a density-matching 
solvent, Fig. [T] We show that the transition from col¬ 
loidal to granular behaviour is driven by the widespread 
phenomenon of shear thickening EH at a critical ‘onset 
stress’, whose scaling with particle size (oc d~ 2 ), Fig. [2] 
differs from that of the intrinsic stress scale (oc d 3 ). Un¬ 
derstanding this complexity leads to a unified description 
of hard particle rheology over the whole range of d, Fig. [3] 
Interestingly, our results also confirm recent theory and 
simulations Ennuis] ascribing shear thickening to the 
formation of frictional contacts. 

We used polymethylmethacrylate (PMMA) parti¬ 
cles stabilised by 5-10 nm poly-12-hydroxystearic acid 
(PHSA) ‘hairs’ QU in a density-matching mixture of 
cyclolreptylbromide (CHB) or cyclohexylbromide and de- 
calin. We present data for d = 404 nm (from static light 
scattering) and 3770 nm (from microscopy) particles with 
polydispersity ~10% (from light scattering and electron 
microscopy, respectively). Data for d 280 nm (in de- 
calin), 268 nm, 912 nm, 1800 nm and 4500 nm (in CHB 
and decalin) give the same picture. Samples were pre¬ 
pared by diluting a close packed sediment, using simu¬ 
lations m to estimate 4>bcp of polydisperse HS. The 
solvent viscosities were r/f = 2.83 rnPa.s for the large and 
2.4 rnPa.s for the small particles at 19°C. Adding an ex¬ 
cess of screening salt tetrabutylammonium chloride did 
not change the rheology; we present salt-free data. 

Rheology was performed in an Anton Paar 301 rheome- 
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FIG. 1. Rheology in the colloidal, intermediate and granular size regimes, (a) Relative viscosity 17 versus shear stress in units 
of Pa and ksT/d 3 for d = 404 nm spheres. At a < 10 Pa we performed downward sweeps in 7 and at er > 10 Pa, upward sweeps 
in <t, at volume fractions (%) as labelled. Solid lines: fits based on [T2]; finely-dotted lines: schematics based on literature data 
(with sparsely-dotted = unstable states). Samples shear thicken above a ^-independent onset stress cr* (vertical dashed line). 
Colour scheme: blue = frictionless interactions (a < a*), black = shear thickening, red = frictional interactions (cr » a*). The 
unshaded region is accessible using our rheometer, which reaches maximum and minimum shear rates of 8000 s -1 and 10~ 3 s -1 
respectively, and a minimum stress of 1CT 2 Pa; the maximum accessible stress is set by a d-dependent fracture stress, crh (b) 
Main: p(</>) for the the limiting high-shear viscosity, 171 , in (a) (■), and the lower, i 7 \((j>) (■), and upper, r/ 2 ((j)) (□), branches in 
(c). Solid red line: least squares fit to 72 (<p) = A(1 - 4>l<j) rn )~ n with A = 0.20(9), = 0.558(5) and n = 2.2(2). 771 (0) data for 

other sizes of PMMA spheres in this work: d = 912 nm (v) and 1800 nm (▼). Other symbols: literature high shear viscosities 
(with <f> shifted by up to 5%) for sterically stabilised PMMA [jS] E; 10 Eli sterically stabilised silica [ 5 ] and glass beads m 
(•) and (o) = lower and upper branches from [17]. Inset: 172 ((f)) versus (0 m ~(j>) including the upper branch from [17] ■ (c) 17 (a) 
for d = 3770 nm spheres, all from upward cr sweeps apart from <j) = 58%, which is from a cr downward sweep. The flow in both 
(a) and (c) was unsteady for (f> > 0.56, and points represent temporal averages, (d) Schematic for the rheology of d = 350 /cm 
spheres calculated from theoretical fits to m in (c). 


ter in truncated cone and plate geometry (cone angle 
1°, radius 25 mm, truncation gap 100 pm) in stress- 
controlled mode unless otherwise stated. A sandblasted 
cone (surface roughness ~ 10 pm) and a base plate rough¬ 
ened with silicon carbide powder (surface roughness ~ 
5 pm) were used. Even with a solvent trap, artefacts 
due to drying were evident if samples were left to equi¬ 
librate at shear rates 7 < 10 “ 2 s -1 . Thus, we worked at 
7 > 0.01 s -1 . 

To establish a baseline, we first explore d = 404 nm col¬ 
loids. Figure[l|i shows the relative viscosity 17 = (a/j)/r)f 


(solvent viscosity 77 / = 2.4 mPa.s) as a function of applied 
shear stress cr. The region in each 17 (a) plot in Fig. [l] not 
shaded dark grey is the ‘observable window’ accessible by 
conventional rheology, and points are experimental data. 
Additionally, the light grey region in Fig. E* is accessible 
but typically not probed in previous work 00 . 

Our data for </> = 0.40 and 0.47 show classic shear thin¬ 
ning to a high-shear Newtonian plateau, rji((j)), in the ob¬ 
servable window. The <fi = 0.53 sample behaves similarly 
in this window, but shear thickens beyond it. At <j> = 0.56, 
we see the onset of discontinuous shear thickening (DST), 
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FIG. 2. Scaling of the onset stress with particle size. On¬ 
set stress a* versus particle diameter d for PHSA-stabilised 
PMMA spheres in this and previous [9] [16] 1T9H22] work. The 
solid line is a least-squares fit to a power law a* = Bd a , with 
B = 2 X nr 11 J.nT 1 1 and a = -1.9(1). 

i.e. the gradient dlogr// log cr reaches 1 , corresponding to 
a vertical flow curve, dj/da = 0. Our ? 7 i(</>) data, Fig.[l]D 
(■), agree with previous measurements of the ‘high-shear 
viscosity’ of various colloids [310161 OH [[7]. Within ex¬ 
perimental uncertainty, shear thickening begins at a fixed 
‘onset stress’, cr* ^ 200 Pa, for all (j> . At a higher </>- 
independent stress erf ^ 10 5 Pa, samples fracture - this 
is the high-cr boundary to the observable window. (The 
lower right, slanted, boundary is due to an inertial insta¬ 
bility expelling samples from the instrument.) 

Increasing d by an order of magnitude to 3770 nm has 
a dramatic effect on the rheology in the observable win¬ 
dow, Fig. [l]:. Now, at 0.41 < </> < 0.54, there is a transi¬ 
tion from a lower to a higher, shear-thickened, Newtonian 
plateau as a increases HDUBJ. The transition to a shear- 
thickened state again begins at a (^-independent onset 
stress, now at cr* ^ 2 Pa. Approximately the same onset 
stress applies at 4> > 0.56, but now the high-shear plateau 
vanishes. The data become noisy and dlogi]/dloga ^ 1, 
which is the DST limit. The onset of fracturing in the 
sample is now at erf ^ 500 Pa, above which the data show 
poor reproducibility and strong history dependence. 

The two plateau viscosities for the d = 3770 nm parti¬ 
cles fall on separate diverging branches, Fig. [TJd (■ and 
□) [23]. The lower plateau viscosity here corresponds 
to 7 /i(</>) for the d = 404 nm colloids (■): both diverge 
at <p ^ 0.64 = (f>R cp- Thus, we infer an experimentally- 
inaccessible shear thinning regime for these particles at 
lower cr, sketched schematically in Fig. [lj:. The high vis¬ 
cosity branch, 772 (/>), diverges at (f> m ^ 0.55 < (j> rcp- In¬ 
terestingly, this is close to random loose packing for fric¬ 
tional granular spheres 124] . We infer a corresponding 
regime of shear-thickening to a plateau viscosity in the 
404 nm particles outside the observable window, Fig. UK 


Two viscosity branches can be extracted from previous 
work on PMMA spheres d; each extracted branch can 
be collapsed onto our corresponding branch by a small 
shift in tp (Fig. m Note that it would be meaning¬ 
less (though mathematically possible) to collapse the two 
branches in Fig. [TJd into a single curve by a large (~ 0.1) 
shift in (j): these branches are observed together in the 
same experiment, and correspond to distinct phenomena. 

Recent theory |T2] and simulations ei Ha predict just 
such a two-branch structure for the viscosity in non- 
Brownian (granular) suspensions, and shear thickening 
is associated with transition from the low to the high- 
viscosity branch above an onset stress cr* (see also [25]). 
At cr « cr*, particles do not touch and all contacts in 
the system are lubricated, while for cr» cr*, all particles 
are pressed into frictional contact. The point at which 
the high-viscosity branch diverges, (/> m , decreases with 
increasing (static) friction coefficient, /t p , between par¬ 
ticles El an( I 0m = 0RCP only if /i p = 0. Our data for 
particles in the transitional size regime, Fig. [l]c, are con¬ 
sistent with this scenario; indeed, the solid lines are fits 
of our data to theory [32] (see Supplemental Material [26] 
for details). 

The onset stress in this theory, cr*, arises physically 
from the presence of barriers, e.g. from PHSA ‘hairs’ 
in our case, that stabilise particles against entering the 
primary van der Waals (vdW) minimum of their mutual 
interaction potential, and therefore prevent interparticle 
contact. When the applied stress exceeds a critical value, 
cr*, however, these barriers are overcome and particles are 
pressed together to make frictional contacts. 

We find, Fig. [2] that cr* decreases with particle size as 
cr* _ Bd a with B = 2 x KT 11 J.m -1 ’ 1 and a = -1.9(1). 
For charged-stabilised colloids, one expects cr* oc d~ 2 [2J; 
but charge is irrelevant in our case, because cr* is un¬ 
changed by adding salt (data not shown). More rele¬ 
vantly, cr* oc d~ 2 [27] and oc d~ 175 jT5| scalings are pre¬ 
dicted theoretically for stabilising polymer brushes, and 
d ~ 2 scaling is found in other sterically-stabilised PMMA 
particles [191. 

If we take our data as supporting cr* oc dr 2 , then a 
constant force /* is needed to push particles into fric¬ 
tional contact: cr* 0 : /*/<j 2 with /* = 3 fcsT/nm. This is 
comparable to the 6 ksT /nm measured between similar 
PHSA hairs in a different geometry [28]. The microscopic 
origins of this constant force /* oc d° are at present un¬ 
clear (see jH] and Supplemental Material [26]), although 
cr* oc d ~ 2 may be generic: a d~ 2 scaling of a critical onset 
shear rate is reported in a review of diverse systems [29] 
(where the substantial data scatter may come from using 
7 rather than cr as the scaling variable). 

The interplay between the d~ 2 scaling of the onset 
stress for shear thickening, cr*, and the d~ 3 scaling of 
the intrinsic stress scale controls the colloidal to granular 
crossover in our system. The progression of observable 
rheology is set out schematically in Fig.[3]and its caption. 
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FIG. 3. Hard-sphere rheology unified. Relative viscosity, 
77 , versus dimensionless stress, a. The observable window 
(shaded) for colloids shows shear thinning, which begins at 
<7 - 1; shear thickening, beginning at a dimensionless onset 
stress of <7* = a*d 3 /ksT = /* d/ksT, occurs towards the high 
end of accessible stresses. As the particle diameter d increases, 
<7* shifts right as d, and the second Newtonian regime in the 
77 ( 17 ) curve is stretched out. However, the (shaded) observ¬ 
able window shifts right faster, as d 3 . Thus, shear thinning 
becomes unobservable for intermediate particles sizes, and the 
shear-thickened state fills the observable window of the largest 
particles (grains). Figure^, c, d show three snapshots of this 
scenario. 


(An alternative summary of our findings in the form of 
a ‘rheological state diagram’ is shown in Supplemental 
Material Fig. SI [2B].) This scenario should be valid for 
any system in which a* oc with A > 0, with the exact 
value of A controlling the sharpness of the crossover. Ex¬ 
cept for A = 3, the different scaling of cr* and the intrinsic 
stress (always oc dr 3 for dimensional reasons) means that 
a single set of master curves cannot be found to describe 
the suspension rheology for all d. 

The existence of an ‘onset stress’ <7* means that resid¬ 
ual vdW attraction sets a practical limit to the largest 
particles we could in principle study. For d > 20 p m, the 
vdW attraction at the point when the PHSA hairs on 
neighbouring particles just touch exceeds a* (see Supple¬ 
mental Material [25] for details), and masks the underly¬ 
ing shear thickening HU: the system is now a particulate 
gel. Thus, our system can reach the threshold but cannot 
properly probe the granular regime of repulsive particle 
rheology. 

While our work is not focussed on shear thickening per 
se, our data impact on the understanding of this ubiqui¬ 
tous phenomenon. One theory, e.g. [3D], ascribes shear 
thickening to hydrodynamic interactions (His) alone, 
with interparticle friction playing no role. Simulating 
a system of frictionless spheres with only His gives a vis¬ 
cosity increase that is weak and continuous m- It is un¬ 
clear how the discontinuous shear thickening we observe 
could arise in this framework, whereas a recent theory 
of friction-driven shear thickening [12] can quantitatively 
describe our results. We also note that the shear thick¬ 
ening we observe is distinct from recent work on inertial 
systems [321 S3] , as our particle Reynolds number is at 


most « 10 2 , and typically < 10 4 at the onset of thick¬ 
ening. 

To conclude, we have shown that the transition from 
colloidal to granular rheology is driven by shear thick¬ 
ening. Our data are consistent with recent suggestions 
El HI [HI that shear thickening is associated with the 
development of frictional particle contacts at an onset 
stress, cr*, which we find to decrease with particle size 
as dT 2 . Thus, particles with d < 1 pm will behave as 
frictionless, Brownian hard spheres at most accessible 
stresses. For intermediate-sized particles, frictionless and 
frictional states are observed at low and high stress re¬ 
spectively. Finally, a particulate suspension is granular 
when cr* is much smaller than commonly-encountered 
stresses; such a suspension is ‘always shear thickened’. 
The size at which this happens depends on the ‘stabilis¬ 
ing force’ /*, and therefore the surface chemistry. 

That the onset stress is readily accessible and routinely 
exceedable for a suspension of particles with 1 pm < 
d < 50 pm, Figs. [I]:' and [2] has significant practical 
consequences, especially for concentrations in the range 
(j> m < <j> < i^rcp- Once cr* is exceeded, there is no fric¬ 
tional branch with finite viscosity to which the suspension 
may jump, Fig. [Tj>. The system shear jams, and shows 
a qualitative change in its rheology consistent with pre¬ 
vious work on concentrated suspensions [34] ■ The flow 
is unsteady, shows edge fracture and wall slip (confirmed 
by imaging), and becomes strongly history-dependent. 

Understanding polydisperse industrial suspensions 
whose particle size distributions span all three regimes re¬ 
mains a formidable challenge. Our unified description of 
particulate rheology over all sizes, summarised in Fig. [3] 
has laid the foundation for this challenge to be met. 
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